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WIND-TUNNEL  RESEARCH  ON  THE  MECHANICS  OF  PLUMES 
IN  THE  ATMOSPHERIC  SURFACE  LAYER 

1.  INTRODUCTION 
1 . 1  Purpose 

The  three-year  research  program  is  intended  to  establish  the 
foundations  for  a  long-range  comprehensive  investigation  of  aero- 
collodial  dispersion  and  evolution  in  the  atmospheric-surface-layer 
(ASL)  environment  through  use  of  simulated  ASL  flows.  It  was  formulated 
to  achieve  the  following  objectives: 


1.  Establishment  of  a  set  of  reference  data  on  mean  flow  and 
turbulence  characteristics  for  the  full-range  of  possible  ASL 
simulations  in  the  Meteorological  Wind  Tunnel  (MWT)  of  the 
Fluid  Dynamics  and  Diffusion  Laboratory  (FDDL) . 

2.  Development  of  dispersion  data  for  the  reference  flows  that 
will  serve  as  reference  plumes  to  evaluate  effects  of  aero- 
colloids  on  dispersion  mechanics. 

3.  Design  and  initiate  studies  to  determine  the  effects  of  aero- 
colloidal  particle  size  and  concentration  on  concentration 
distributions  and  surface  deposition. 

4.  Design  and  develop  a  capability  for  research  on  aerocolloidal 
kinetics  in  the  MWT. 

5.  Establish  the  foundation  for  a  long-range  research  program  on 
plume  mechanics. 


The  proposed  research  objectives  are  further  described  in  specific 
tasks  planned  for  each  of  the  three  years. 

1 . 2  First-Year  Tasks 

Initial  efforts  of  the  research  program  were  directed  to  the 
development  of  reference  data  for  a  wide  range  ASL  simulations  in  the 
MWT.  Mean  flow,  turbulence  characteristics  and  dispersion  data  were 
investigated  and  carefully  documented.  Additionally,  the  design  of 
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second-year  experiments  with  aerocol loida  1  particles  was  initiated. 
These  three  first-year  objectives  are  expanded  upon  in  the  task  des¬ 
criptions  which  follow: 

Task  1 -A :  Establish  characteristics  of  meteorological  variables  for 
simulated  ASL's  in  the  MWT  that  include  full  range  of  the  wind 
tunnel  capability  to  simulate  stable  and  unstable  thermal 
stratification  beginning  at  ground  level  and  elevated  inver¬ 
sions.  Careful  documentation  of  the  mean  flow  and  turbulence 
characteristics  is  required  to  provide  basic  reference  flow 
fields  for  the  plume  dispersion  studies. 

A.  Flow  Variables 

1.  Lower  boundary  --  plane 

2.  Boundary  roughness  --  2  cases  (smooth,  1  rough) 

3.  Thermal  stratification  —  4  cases  (neutral,  1  stable, 

1  unstable,  1  elevated  inversion) 

B.  Measurements  (for  each  combination  of  flow  variables) 

1.  Mean  velocity  and  temperature  —  5  vertical  profiles; 

along  centerline  at  x  =  15,16,17,19, 
22  m  from  test-section  entrance 

2.  Turbulence  intensities  --  5  vertical  profiles;  same 

location  as  for  mean  velocity  and 
temperature _  _  _ 

3.  Velocity-velocity  correlations  (u'w' ,  u'v' ,  v'w' )  --  1 

vertical  profile,  x  =  16  m  on 
centerline 

4.  Velocity  power  spectra  (u* ,  v' ,  w')  --  5  points;  x  =  16  m 

5.  Turbulence  scales  (longitudinal,  lateral,  vertical)  --  5 

points;  x  =  16  m 

Task  1-B:  Develop  a  set  of  reference  plume  characteristics  for  the 
reference  flow  fields.  The  reference  plumes  are  intende  1 
to  provide  basic  characterization  of  turbulent  diffusion 
by  atmospheric  turbulence  without  modifications  caused  by 
the  diffusing  material.  Reference  plumes  developed  by  a 
single  source  emission  at  several  distinct  elevations  in 
the  simulated  ASL,  including  ground  level,  are  proposed 
for  reference  documentation. 
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A.  Flow  Variables:  same  as  Task  1-A 

B.  Materi.il  Source  Variables 

1  .  IT  me :  s  t  eady 

2.  Spatial:  point;  3  elevations  at  x  =  16  m 

3.  Material:  noil-reacting  gas,  with  hydrocarbon  tracer, 

neutral  buoyancy  relative  to  atmospheric 
density  at  release  elevation 

A.  Momentum:  near  zero 

C.  Measurements 

1.  Mean  concentration  of  source  material:  distributions  to 
describe  cross  section  of  plume  at  4  locations  downwind  of 
source 

Task  1-C:  Design  experiments  to  investigate  plume  mechanics  for 
stable  aerocol loidal  particles  in  the  simulated  ASL. 

A  plume  study  to  be  performed  in  Year  2  will  be  designed  to  include 

the  following  general  specifications  -- 

A.  Flow  Variables:  same  as  Task  1-A 

B.  Material  Source  Variables: 

1.  Time:  steady 

2 .  Spat  ia  1  :  point 

3.  Material  (tentative):  monodisperse  polystyrene  latex 

spheres  (0.5  -  5.0  pm) 

C.  Measurements 

1.  Mean  concentration  distributions  in  plume 

2.  Surface  deposition 

1 . 3  Report  Organizat ion 

The  MWT  was  operated  in  several  different  configurations  to  achieve 
the  thermal  stratifications  and  boundary  roughnesses  specified  in  the 
preceding  Tasks  1-A  and  t-B.  These  configurations,  as  well  as  descrip¬ 
tions  of  the  test  instrumentation  and  measurement  apparatus,  are  re¬ 
corded  in  the  next  section. 

Experiments  were  conducted  to  obtain  reference  data  for  three 
differing  atmospheric  stabilities.  The  results  of  these  experiments  are 
discussed  in  Sections  3.0,  4.0  and  5.0  for  neutral,  stable  and  unstable 
flow,  respectively.  Each  section  adheres  to  a  similar  format.  The 
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velocity  profiles  obtained  above  a  smooth  floor  are  first  described, 
followed  by  the  velocity  fields  recorded  over  the  rough  floor.  The 
final  part  of  each  section  contains  a  discussion  of  the  measured  con¬ 
centration  fields  for  both  roughness  conditions.  Tabulated  data  and 
graphic  comparisons  used  in  the  presentations  are  contained  in  tables 
and  figures  located  at  the  end  of  each  section. 

Section  6.0  contains  Conclusions  and  Commentary  from  the  previously 
discussed  experiments. 

The  appendices  to  this  report  contain  development  of  special  topics 
entitled  "Wind-Tunnel  Simuluation  of  Particle-Plumes  in  the  Atmospheric 
Surface  Layer"  and  "Length-Scales  Describing  the  Vertical  Diffusion  from 
Ground-Level  Sources". 
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2.  EXPERIMENTAL  PROCEDURES 
2 . 1  Wind-Tunnel  Conf igurat ion 

All  experiments  were  conducted  in  the  MWT  at  CSU.  Design  and 
operation  of  the  wind  tunnel  are  described  in  detail  by  Cermak  (1981). 
Elevation  and  plan  views  of  the  MWT  are  shown  in  Figure  2-1. 

Figure  2-2  is  a  three-dimensional  representation  of  the  entrance  to 
the  MWT  test  section.  All  four  surfaces  immediately  preceding  the  test 
section  were  lined  with  gravel.  The  sawtooth  shaped  trip  which  framed 
the  entrance  was  also  used  as  a  reference  for  locating  the  point  source 
within  the  test  section. 

Four  evenly  spaced  91  cm  tall  spires  were  placed  across  the  tunnel, 
next  to  the  sawtooth  trip,  to  simulate  the  desired  ABL  within  the  test 
section.  Pertinent  theories  of  ABL  simulation  are  discussed  in  detail 
by  Cermak  (1971,  1982).  Figure  2-3  contains  a  close-up  pictorial 
representation  of  the  MWT  entrance  configuration,  while  Figure  2-4 
contains  a  dwonwind  view  of  the  entire  MWT  test  section,  without 
roughness . 

A  portion  of  the  twelve  auxiliary  cooling  plates,  which  covered  one 
third  of  the  floor  upwind  of  the  test  section,  is  visible  in  the  photo 
reproduced  in  Figure  2-5.  These  plates,  which  remained  in  place  during 
tests  in  all  three  stratifications,  have  a  slightly  rougher  surface  than 
the  permanently  installed  aluminum  plates,  which  make  up  the  next 
12.28  m  of  tunnel  floor. 

A  change  in  the  aerodynamic  surface  roughness,  zq,  from  a  lower 
value  for  the  'smooth'  case  to  a  higher  value  for  the  'rough'  case  was 
achieved  by  the  installation  of  1.6  cm  twisted-link  chain  elements 
placed  21.6  cm  apart  along  the  entire  test  section  of  the  MWT.  The 
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origin  of  the  coordinate  system  used  to  reference  the  velocity  and 
concentration  measurements  was  established  at  the  location  of  the 


ground-level  point  source  (16.12  m  downwind  of  the  sawtooth  trip). 
Figure  2-6  shows  the  arrangement  of  the  roughness  chains. 

A  schematic  drawing  summarizing  the  configuration  of  the  MWT  is 
presented  in  Figure  2-7.  This  drawing  also  shows  the  location  of  the 
measurement  stations  which  are  described  in  the  next  sub-section. 

2.2  Test  Measurements  and  Instrumentation 


The  neutral  boundary  layer  was  simulated  in  the  MWT  with  existing 
ambient  conditions  whereas  stratified  boundary  layers  were  simulated  by 
controlling  the  environment  in  the  wind-tunnel.  Table  2-1  presents  the 
mean  velocity  and  temperature  conditions  maintained  in  the  MWT  during 
experiments  with  each  of  the  three  stabilities.  Free-stream  velocity 
and  temperature  were  continuously  monitored  with  a  pitot  tube  and 
thermistor  positioned  1.09  m  above  the  floor,  and  0.85  m  upwind  of  the 
source  (see  Figures  2-6  and  2-7). 

To  create  the  stable  boundary  layer,  an  ethylene-glycol  solution 
was  circulated  through  the  auxiliary  cooling  plates  and  the  permanent 
aluminum  plates  to  chill  the  floor,  while  the  air  circulating  in  the 
closed-loop  tunnel  was  heated. 

The  unstable  boundary  layer  was  simulated  by  electrically  heating 
the  aluminum  plates  and  simultaneously,  cooling  the  air  circulating  in 
the  MWT . 

Using  the  location  of  the  ground-level  source  as  a  coordinate 
reference,  velocity,  turbulence  and  temperature  data  for  each  stability 
and  surface  roughness  condition  were  measured  at  downwind  distances  of 


0,  2.80  and  6.05  m  along  the  centerline.  These  measurements  were 


obtained  with  a  constant  temperature  cross-film  anemometer,  TS1  model 
1241-10,  ami  a  iesi statue  wire  temperature  sensor  mounted  in  close 
proximity  to  each  other.  Prior  to  taking  data,  the  sensors  were  cali¬ 
brated  in  a  temperature  and  velocity  controlled  enclosure.  Empirical 
equations  were  developed  to  relate  anemometer  voltage,  velocity  and 
temperature.  Using  these  equations,  an  on-line  computer  permitted 
nearly  instantaneous  output  of  the  velocity  and  temperature  data 
measured  in  the  MWT . 

Two  orientations  of  the  cross-film  furnished  velocity  and  turbulence 
data  along  all  three  primary  axes.  Mean  and  cross-correlation  values 
were  computed  from  600-1000  instantaneous  readings  per  second,  recorded 
over  a  period  of  60-120  seconds  with  sampling  rate  and  duration 
dependent  upon  the  relative  stability  of  test  conditions. 

Neutrally  buoyant  hydrocarbon  tracers  were  released  from  0.635  cm 
diameter  brass  tubes  (referred  to  elsewhere,  and  here,  as  the  point 
sources),  located  at  heights  of  zS  =  0,  10.00  and  20.0  cm  above  the 
floor.  The  effective  height  of  the  ground-level  source  was  one-half  the 
tube  diameter,  or  approximately  3  mm.  Figure  2-8  illustrates  the 
arrangement  of  the  point  sources. 

The  tracers  (9.98%  ethane,  4%  CO^,  86.0%  or  5.6%  propane,  94.38% 
N^)  were  released  through  the  sources  with  an  exit  velocity  ~  88  cm/sec 
while  being  maintained  at  the  average  temperature  existing  at  source 
he i ght . 

Tracer  com ent rat  ions  were  measured  at  distances  of  0.85,  1.95, 
2.80,  3.90  and  6.05  m  downwind  from  the  source  for  all  three  stabil¬ 
ities.  Rapid  diffusion  in  the  unstable  flow  dictated  additional 
measurements  at  0.45  and  1.35  m. 


Samples  ol  the  tracer  were  drawn  by  a  near- isokinetic  sampling 
system  through  1.59  mm  I.D.  tubes  into  a  gas  chromatograph  analyzer. 
The  sampling  tubes,  mounted  on  a  rake,  as  shown  in  Figure  2-9,  provided 
a  grid  for  measuring  diffusion  within  the  tracer  plume.  Up  to  40  five- 
minute  tracer  samples  were  simultaneously  collected  by  the  sampler, 
shown  in  Figure  2-10a,  for  analysis  int  he  Hewlett-Packard  gas 
chromatograph  (GC),  shown  in  Figure  2-10b. 

The  tracer  gas  data  obtained  from  the  GC  analyses  wer  normalized 
using  the  equation 


N.C.  = 


(2-1) 


where  N.C.  =  normalized  concentration  (m  ) 

X  =  measured  concentration  (ppm) 

u^  =  free-stream  velocity  (m/s) 

Xq  =  source  strength  (ppm) 
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V  =  volume  flowrate  (m  /s). 

All  diffusion  data  included  in  this  report  are  in  this  normalized 


2 . 3  Flow  Visualization 

Titanium  dioxide  (TiO^)  'smoke'  was  released  from  the  point  sources 
for  each  stratification  to  observe  plume  behavior  and  to  ascertain  the 
normal  operation  of  lthe  MWT.  Plume  diffusions  for  the  three  stratifi¬ 
cations  and  three  source  heights  were  documented  on  16  mm  movie  film  and 
35  mm  photos.  Representative  plumes  from  a  ground-level  source  for 
stable,  neutral  and  unstable  atmospheric  conditions  are  presented  in 
F i gu res  2-11  and  2-12. 


Table  2-1.  Temperature  and  Velocity  Conditions  in  the  MWT  for  Plume 
Mechanics  Experiments 


Atmospheric 

Condition  Roughness 


smooth 

23.2 

3.00 

(22.5  to  24.0) 

(2.88  to  3.06) 

Neutral 

rough 

- 

23.5 

2.99 

(22.5  to  25.0) 

(2.91  to  3.05) 

smooth 

~±1 

77.0 

2.45 

(75.7  to  79.3) 

(2.43  to  2.47) 

Stable 

rough 

~±1 

78.7 

2.44 

(73.0  to  80.3) 

(2.39  to  2.51) 

smooth 

157 

7.7 

2.43 

Unstable 

(7.4  to  8.4) 

(2.40  to  2.46) 

rough 

146 

7.5 

2.44 
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Figure  2-5.  Auxiliary  cooling  plates  which  covered  upwind 
one-third  of  MWT  test-section. 


Figure  2-6.  Upwind  view  of  MWT  test-section  with  roughness  in  place 


ire  2-7.  Schematic  of  MWT  configuration  and  measurement  locations 


Figure  2-10.  Photographs  of  (a)  the  Tracer  Sampling  System  and 
(b)  the  HP  Integrator  and  Gas  Chromatograph. 
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NEUTRAL  FLOW 


3.1  Velocity  Field  Over  the  Smooth  F loor 

Mean  velocity  distributions  u(z)  were  measured  at  distances  of 
x  =  O.C,  2.80  and  6.05  m  with  respect  to  the  source.  These  velocity 
profiles  are  presented  in  Figure  3-1  and  indicate  good  horizontal  homo¬ 
geneity.  However,  since  the  flow  transitions  from  a  slightly  rough  to 
smooth  surface,  in  accordance  with  Antonia  and  Luxton's  (1972)  findings, 
the  shear  stress  in  the  central  region  of  the  boundary  layer  increases 
while  the  surface  shear  stress  is  reduced.  With  increasing  distance 
from  the  transition  of  the  surface  roughness,  the  surface  shear  stress 
gradually  increases.  These  effects  can  be  seen  in  Figure  3-2  which 
presents  the  shear  stress  measurements  at  distances  of  x  =  0.0,  2.80  and 
6.05  m  from  the  source. 

The  average  value  of  the  shear  velocity  for  this  flow  in  the  region 

of  interest  is  estimated  from  Figure  3-2  to  be  u*  =  0.09  m/sec.  The 

measured  longitudinal  turbulence  intensity  (u  /U)  and  vertical  turbu- 

rms 

lence  intensity  (Wrms/U)  distributions  are  presented  in  Figures  3-3  and 
3-4  for  measurements  made  at  distances  of  x  =  0.0,  2.80  and  6.05  m  from 
the  source. 

From  the  velocity  and  turbulence  data  for  neutral  flow  over  the 
smooth  floor  it  may  be  concluded  that  the  boundary  layer  developed  in 
the  wind  tunnel  was  not  a  perfect  model  of  an  ideal  neutral  boundary 
layer.  However,  its  major  characteristics  are  typical  or  many  real 
neutral  atmospheric  boundary  layers. 

3 . 2  Velocity  Field  Over  the  Rough  Floor 

The  mean  velocity  distributions  u(z)  measured  at  distances 
of  x  =  0.0,  2.80  and  6.05  m  with  respect  to  the  source  are  presented  in 


33 


Figure  3-5.  These  velocity  profiles  also  indicate  good  horizontal 
homogeneity.  The  transition  in  surface  roughness  in  this  case  is 
comparatively  smaller  than  that  for  the  case  of  flow  over  the  smooth 
surface.  Hence  distributions  of  the  shear  stress  at  distances  of 
x  =  0.0,  2.80  and  6.05  m  from  the  source,  presented  in  Figure  3-6 
exhibit  the  familiar  monotonica 1 ly  decreasing  behavior  throughout  the 
central  region  of  the  boundary  layer. 

The  average  shear  velocity  for  this  flow  in  the  region  of  interest 
is  estimated  from  Figure  3-6  to  be  u*  =  0.135  m/sec.  The  longitudinal 
turbulence  intensity  and  vertical  turbulence  intensity  distributions  are 
presented  in  Figures  3-7  and  3-8  for  measurements  made  at  distances  of 
x  =  0.0,  2.80  and  6.05  m  from  the  source. 

Table  3-1  compares  some  properties  of  the  simulated  wind  flow  in 
the  wind  tunnel  with  Robins'  (1978)  data,  and  with  typical  values  of 
these  properties  found  in  the  atmospheric  boundary  layer.  In  the  table, 
zq  denotes  the  roughness  length,  u*  the  friction  velocity,  Ure£  the 
reference  velocity,  u,  v  and  w  the  turbulence  components  in  the  longi¬ 
tudinal,  lateral  and  vertical  directions,  and  n  the  exponent  of  the 
velocity  profile: 
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the  surface.  The  measurements  were  made  at  downstream  distances  of 
x  =  0.85,  1.95,  2.80,  3.90  and  6.05  m  from  the  source.  The  data 

corresponding  to  the  ground-level  source  were  fitted  to  the  empirical 
semi-equation 


(3-2) 


q — =  exp[  (-«n2)(^^)2]exp[  (-£n2)(^-)S]  , 

mo  v  *  ‘ 


where  the  knowns  are  C:  normalized  observed  concentration,  m 

y,z:  measuring  location,  cm 
and  the  unknowns  are 

_2 

C  (x):  maximum  normalized  plume  concentration,  m 
max  ’ 

dy:  horizontal  drift  of  maximum  concentration,  cm 
Ay:  horizontal  plume  spread  constant,  cm 
Az:  vertical  plume  spread  constant,  cm 
s:  vertical  plume  spread  exponent. 

Estimates  of  the  five  unknowns  in  the  above  equation  for  the  five 
downstream  measurement  stations  are  presented  in  Table  3-2. 

Data  collected  in  the  horizontal  and  vertical  planes  corresponding 
to  the  source  located  at  a  height  of  10  cm  above  the  surface  were  fitted 
to  the  semi-empirical  equations 


“So  = 


£ — =  exp((-£n2)(“)2]+exp[(-£n2)(~)2]  , 
max^ 


(3-3) 


(3-4) 


where  H  is  the  known  source  height  above  the  surface.  Estimates  of  the 

four  unknowns  C  (x),  dv.  Ay  and  Az  in  the  above  equations  for  the  five 
max 

downstream  measurement  stations  are  presented  in  Table  3-3. 


Equations  for  diffusion  in  the  horizontal  plane  from  the  ground 
level  as  well  as  the  elevated  sources  were  based  on  a  Gaussian  diffusion 
model.  Figures  3-9  and  3-10  show  typical  horizontal  concentration 
profiles,  measured  as  well  as  predicted  by  these  equations  using  the 
computed  values  of  the  unknowns,  downstream  from  the  ground-level  and 
elevated  sources  respectively.  These  figure  clearly  demonstrate  that 
the  Gaussian  model  is  suitable  for  horizontal  diffusion. 

Diffusion  in  the  vertical  plane  is  fairly  well  described  by  the 
Gaussian  model  for  the  case  of  the  elevated  source,  as  can  be  seen  in 
Figure  3-11,  but  this  model  fails  to  accurately  describe  the  diffusion 
in  the  vertical  plane  from  a  ground-level  source.  Hence  for  the  ground- 
level  source,  diffusion  in  the  vertical  plane  is  described  by  a  non- 
Gaussian,  variable-exponent  model.  Figure  3-12  illustrates  the 
suitability  of  an  exponential  model  in  the  present  investigation. 
Non-Gaussian  distributions  of  tracer  material  were  also  observed  in  the 
diffusion  experiments  of  the  Prairie  Grass  Project  (Barad  and  Haugen, 
1959)  and  the  Green  Glow  program  (Barad  and  Fuquay,  1962). 

Diffusion  data  in  the  vertical  plane  from  the  elevated  source,  for 
flow  over  both  surface  roughnesses  considered,  shows  that  near  the 
source  the  maximum  concentration  occurs  at  the  height  of  the  source, 
then  it  rapidly  approaches  the  ground  with  increasing  distance  from  the 
source. 

Figures  3-13,  3-14,  3-15  and  3-16  show  the  distributions  of  \^(x), 
Oy(x),  Az(x),  and  <Jz(x)  for  the  two  source  heights  and  two  surface 
roughnesses  considered.  The  o^(x)  and  az(x)  are  standard  deviations  of 
the  horizontal  and  vertical  plume  spreads,  respectively.  The  variation 
of  the  normalized  peak  concentrations,  with  distance  downstream  from  the 
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Table  3-1.  Comparison  of  Simulated  and  Atmospheric  Flows 


Category 

6 

(m) 

2 

o 

(m  at 
full 
scale) 

n 

u  */U  , 

ref 

u  / 1 

rms 

Present 

data 

Smooth 

0.95 

0.01 

0.15 

0.037 

2.1 

Rough 

0.90 

0. 16 

0.25 

0.055 

2.5 

Robins ' 

data 

Rural 

0.6 

0.03 

0.14 

0.045 

2.4 

Urban 

2.0 

1.3 

0.22 

0.063 

2.6 

Atmosphere# 


Rura  1 

600 

0.01-0.07 

1/7 

0.037-0.045 

2.5 

Suburban 

600 

1.0-1. 5 

1/5 

0.045-0.055 

2.5 

Urban 

600 

1. 0-4.0 

1/4 

0.040-0.060 

2.5 

//Source:  Counihan  (1975) 


Smooth 


Distance  Downstream  from  Source  (m) 


(cm)  5.54  20.10  21.07  21.64  28.89  Smooth 
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Figure  3-1.  Mean  velocity  profiles,  smooth  floor. 
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figure  3-10.  Horizontal  concentration  profiles  for  a  10  cm  source 
in  neutral  flow. 
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3-11.  Vertical  concentration  profiles  for  a  10 
in  neutral  flow. 
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Figure  3-12.  Vertical  concentration  profiles  for  a  ground-level  source 
in  neutral  flow. 
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Figure  3-18.  Normalized  peak  concentrations  for  two  source  heights,  smooth  floor. 
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STABLE  FLOW 


4 . 1  Velocity  Field  over  the  Smoo t h  Floor 

Temperature  profiles  T(z)  were  measured  at  distances  of  x  =  0.0, 
2.80  and  6.05  m  from  the  source.  These  profiles  are  presented  in 
Figure  4-1  and  show  good  horizontal  homogeneity. 

Mean  velocity  distributions,  u(z),  measured  at  distances  of 

x  =  0.0,  2.80  and  6.05  m  from  the  source  are  presented  in  Figure  4-2. 

These  profiles  also  show  good  horizontal  homogeneity,  however,  they  also 
indicate  that  the  flow  was  not  ideal.  Such  a  behavior  of  the  flow 
field,  which  is  typical  of  real  atmospheric  flows,  can  be  attributed  to 
the  step  change  in  surface  roughness  from  rough  to  smooth,  as  discussed 
in  an  earlier  section  on  neutral  flows.  Under  these  conditions,  the 
selection  of  idealized  parameters  such  as  the  shear  velocity,  u* ,  the 
surface  roughness,  zo>  and  the  Monin-Obukhov  length  scale,  L,  is 

affected . 

Shear  measurements  made  at  distances  of  x  =  0.0,  2.80  and  6.05  m 
from  the  source  are  presented  in  Figure  4-3.  These  profiles,  too, 
clearly  show  the  effect  of  flow  transition  from  rough  to  smooth  surface. 
Since  stable  flows  are  very  strongly  affected  by  the  transition  in 

surface  roughness,  the  average  value  of  the  shear  velocity,  uvV ,  for  this 
flow  in  the  region  of  interest  is  estimated  by  fitting  data  near  the 
surface  to  the  log-linear  equation, 

u(z)  =  jp  |£n(  — )]  ,  (4-1) 

N  Z 

O 

where  u(z)  is  the  mean  velocity  at  elevation  z,  k  is  von  Karman's 
constant,  and  z^  is  the  aerodynamic  surface  roughness.  The  value  of  the 


shear  velocity  thus  obtained  is  u-  =0.11  m/s. 


From  the  similarity  theory  of  Monin  and  Obukhov,  which  is  based  on 
the  assumption  that  the  flow  is  plane  homogeneous  and  that  vertical 
fluxes  are  constant,  the  following  velocity,  temperature  and  length 
scales  are  obtained: 

The  shear  velocity  u*: 

u"  =  <VPo)!*  (4-2) 

The  friction  temperature  T*: 


T*  =  -H  / (p  C  Ku*) 

o  o  p 


(4-3) 


The  stability  length  L: 


H 


L  =  -U*3/(K  I”  > 

a  o  p 


(4-4) 


where  p  is  the  density,  t  is  the  wall  shear  stress,  H  is  the  wall 

ri  J  r\  7 


heat  flux,  and  g/T  is  the  stability  parameter. 

d 


For  stable  flow  over  a  smooth  surface,  the  friction  temperature  is 
estimated  from  Figure  4-1  to  be  T*  =  11C°.  Combining  equations  4-3  and 
4-4,  the  stability  length  is  expressed  as 


L  =  u*2/(K2TS-  T*)  .  (4-5) 

a 

Using  the  above  equation,  L  is  approximately  16  cm  for  this  flow. 
Figure  4-4  shows  the  variation  of  the  point  Richardson  number  with  z/L 
in  this  case. 

Distributions  of  the  longitudinal  turbulence  intensity  (u  /U)  and 

rms 

the  vertical  turbulence  intensity  (w  /U) ,  measured  at  distances  of  x  = 

J  rms  ’ 

0.0,  2.80  and  6.05  m  from  the  source,  are  presented  in  Figures  4-5  and 
4-6. 


4 . 2  Velocity  Field  over  the  Rough  Floor 

Temperature  profiles  T(z)  measured  at  distances  of  x  =  0.0,  2.80 


and  6.05  m  from  the  source  are  presented  in  Figure  4-7. 


Mean  velocity  distributions,  u(z),  measured  at  distances  of  x  =  0.0, 


2.80  ami  6.05  m  from  tlie  source  are  presented  in  Figure  4-8.  These  pro¬ 
files  show  good  horizontal  homogeneity  and  resemble  profiles  for  ideal 
Hows.  As  discussed  in  the  section  on  neutral  flow  over  a  rough  sur¬ 
face,  in  this  flow,  too,  the  change  in  surface  roughness  near  the  region 
of  interest  is  considerably  smaller  in  comparison  to  the  case  of  flow 
over  the  smooth  surface.  Hence  the  shear  stress  profiles  presented  in 
Figure  4-9  for  measurements  made  at  distances  of  x  =  0.0,  2.80  and 
6.05  m  from  the  source  more  closely  resemble  profiles  obtained  in  ideal 
flows.  The  average  shear  velocity  for  this  flow  in  the  region  of 
interest  is  estimated  to  be  u*  =  0.13  m/s  from  curve  fitting  data  to 
equation  4- 1 . 

The  friction  temperature  for  this  flow  is  estimated  from  Figure  4-7 
to  be  T*  =  14C°.  The  stability  length  is  calculated  from  equation  4-5 
to  be  L  ~  25  cm.  Figure  4-10  shows  the  variation  of  the  point 
Richardson  number  with  z/L  for  this  flow. 

The  longitudinal  and  vertical  turbulence  intensity  distributions 
are  presented  in  Figures  4-11  and  4-12,  respectively,  for  measurements 
made  at  distances  of  x  -  0.0,  2.80  and  6.05  m  from  the  source. 

4 . 3  Measured  Concentration  Field 

In  stable  flow,  diffusion  data  were  measured  for  two  surface 
roughness  conditions,  for  a  continuous  point  source  located  at  eleva¬ 
tions  of  ~0 . 0 ,  10.0  and  20.0  cm  above  the  floor  of  the  wind  tunnel. 
Concentration  measurements  were  made  at  distances  ol  x  -  0.85,  1.95, 
2.80,  3.90  and  6.05  m  downstream  of  the  source  location.  Data  corres¬ 
ponding  to  the  ground-level  source  were  fitted  to  the  semiempi ri t a  1 

equation  3-2  which  yielded  the  five  unknowns  C  tx),  d  ,  \  ,  K  and  s, 
1  J  m.i  x  v  V  / 
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I  I. sit'd  in  Table  4-1,  for  the  various  measurement  distances  from  the 
source . 

Concentration  data  corresponding  to  the  two  elevated  sources  were 

fitted  to  semi emp i r i c a  1  equations  3-3  and  3-4.  Estimates  of  the  four 

unknowns  C  (x),  d  ,  A  and  A  obtained  from  these  regressions,  for  the 
max  y  y  z  & 

five  downstream  measurement  stations,  are  presented  in  Tables  4-2  and 
4-3. 

As  observed  in  the  neutral  case,  Figures  4-13,  4-14  and  4-15 

demonstrate  the  adequacy  of  the  Gaussian  model  to  accurately  describe 
lateral  diffusion  in  the  stable  case,  regardless  of  the  surface  rough¬ 
ness  or  source  height.  Figures  4-16  and  4-17  further  show  the  suit¬ 
ability  of  the  same  model  in  describing  diffusion  (in  the  vertical 
plane)  from  elevated  sources.  Vertical  diffusion  from  the  ground  level 
source  in  stable  flow,  as  for  the  neutral  flow,  is  modelled  by  a  non- 
Gaussian,  variable  exponent  model.  Figure  4-18  illustrates  the 
suitability  of  this  model. 
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Figure  4-1.  Temperature  profiles,  smooth  floor 


MEAN  VELOCITY  CCM/SEO 


Figure  4-3.  Vertical  distribution  of  shear  stress  (u'w')  ,  smooth  floor 
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(•]  gure  4-18.  Velocity  concent  rat  ion  profiles  for  i  ground- level 
Source  in  stable  flow. 


5. 


UNSTABLE  FLOW 


5 . 1  Velocity  Field  over  the  Smooth  Floor 

Measurements  of  the  temperature  profiles  at  the  position  of  the 
source  (x  =  0)  and  at  x  =  2.80  m  and  x  =  6.05  m,  are  shown  in 
Figure  5-1. 

The  measurements  reveal  an  inversion  of  the  temperature  gradient  at 
a  height  h  above  the  ground.  The  height  h  is  slightly  increasing  in  the 
direction  of  the  flow,  from  h  =  0.55  m  at  x  =  0  to  h  =  0.85  m  at  x  = 
6.05  m.  The  value  h  =  0.75  m,  at  x  =  2.80  m  appears  to  be  a  repre¬ 
sentative  height  of  the  unstable  layer  for  the  smooth-floor  configura¬ 
tion.  Using  this  value,  a  specific  heat  flux  w'6'  =  0.63°K-m/sec  and 
an  average  temperature  0  =  283°K,  the  representative  value  of  the  con¬ 
vective  velocity  is  estimated  to  be  approximately  w*  =  0.25  m/sec. 

The  measured  mean  velocity  distributions  u(z)  at  x  =  0,  x  =  2.80  m 
and  x  =  6.05  m  are  presented  in  Figure  5-2.  A  mild  horizontal  non¬ 
homogeneity  is  observed  indicating,  as  shown  by  Antonia  and  Luxton 
(1972),  that  the  boundary  layer  as  a  whole  recovers  very  slowly  after  a 
rough-to-smooth  change  in  the  surface  roughness.  Nevertheless,  both  the 
horizontal  and  the  vertical  nonhomogeneities  in  the  convective  layer  are 
mild  and  the  dimensionless  distance  X  =  xW*/ (Uh)  in  this  layer  may  be 
calculated  with  sufficient  accuracy  using  the  average  value  U 
=  1.90  m/sec.  According  to  Antonia  and  Luxton,  the  shear  stress  after  a 
rough-to-smooth  change  in  the  surface  roughness,  increases  throughout 
the  central  region  of  a  boundary  layer  instead  of  monotonical ly  decreas¬ 
ing.  On  the  other  hand,  the  surface  shear  stress  is  drastically  re¬ 
duced,  then  increases  at  a  lower  rate,  with  x.  The  shear  stress  mea¬ 
surements  at  x  =  0,  2.80  m  and  6.05  m  which  are  presented  in  Figure  5-3 
show  similar  behavior. 


The  average  value  of  the  shear  velocity  in  this  region  is  estimated 
from  figure  5-3  to  he  u-  -  0.085  m/sec,  which  implies  that  the  Monin- 
Obukhov  length  for  this  flow  was  on  the  order  of  0.07m  and  -h/L  ~  10.7. 
The  Mon i n-Ohukhov  length,  L,  which  is  a  measure  of  the  thickness  of  the 
surface  layer  where  the  shear  stress  plays  a  dominant  role,  is  defined 

3  3 

L  =  -(u*)3  0/k  g  wW  ■  (5-1) 

Since  in  this  experiment  the  shear  stresses  in  the  central  region  of  the 
boundary  layer  were  slightly  larger  than  the  shear  stress  at  the  wall, 
it  is  plausible  that  the  surface  layer  extends  beyond  the  above  esti¬ 
mated  height,  or  in  other  words,  the  effective  value  of  -h/L  is  smaller. 

2  ,  2 

The  measured  normalized  vertical  velocity  fluctuations,  w'  /  (w~)  , 
are  plotted  in  Figure  5-4.  The  values  for  the  flow  over  the  smooth 
floor  vary  between  0.1  to  0.3  and  are  of  the  same  order  of  magnitude  as 
the  field  measurements  in  the  Ashchurch  and  Minnesota  experiments  (see 
Lamb,  1981,  Figure  4.18). 

It  may  be  concluded  that  although  the  unstable  boundary  layer 
developed  in  the  wind  tunnel  over  the  smooth  floor  is  not  a  perfect 
model  of  an  ideal  CBL,  its  major  characteristics  are  typical  of  many 
real  convective  atmospheric  boundary  layers. 

5.2  Velocity  F i el d  over  the  Rough  Floor 

The  temperature  distributions  for  this  case  are  shown  in  Figure 
5-5.  The  values  of  h  at  the  different  stations  are  slightly  higher  for 
this  case  and  the  representative  value  of  h  is  estimated  to  he  h  = 


Tin*  mean  velocity  profiles,  shown  in  Figure  5-6,  indicate  that  the 
horizontal  nonhomogene i ty  of  the  velocity  field  in  this  case  is  smaller, 


apparently  because  the  relative  change  in  the  surface  roughness  was 
smaller.  Distribution  of  shear  stresses  at  the  different  stations, 
shown  in  Figure  5-7,  are  also  more  similar.  Using  these  data,  the 
average  horizontal  velocity  for  this  flow  is  estimated  to  be  0  =  1.70 
P’/sec  and  the'  average  shear  velocity  is  estimated  to  be  v*  =  0.115 
m/sec.  Thus,  the  convective  velocity  scale  and  the  Monin-Obukhov  scale 
for  the  flow  over  the  rough  floor  are  of  the  order  of  w*  =  0.28  m/sec 
and  1,  =  0.16  m.  One  may  conclude  that  the  added  roughness  increased  the 
thickness  of  the  surface  layer  and  reduced  the  convective  nature  of  the 
flow.  The  ratio  -h/L  for  this  case  is  smaller  than  six  and  the  vertical 
velocity  fluctuations  have  also  decreased,  as  seen  from  Figure  5-4. 
5.3  Measured  Concentration  Field 

More  than  40  concentration  profiles  were  measured  in  the  study. 
Hepresentat ive  dimensionless  vertical  concentration  profiles, 
Of x ,0 , z)lJh^/Q  plotted  versus  z/h,  are  shown  in  Figures  5-8a  through  5-8£ 
The  concentration  profiles,  Figures  5-8a,b,c,  measured  immediately 
downwind  from  the  ground- level  source  have  an  approximate  Gaussian  shape 


The  profiles  are,  however,  better  described  by  the  function 
C  =  C ( 0 )  exp  [-p(z/A?)S]  , 


(5-2) 


where  fl  =  £n(2)  and  s  is  a  power  smaller  than  2.  The  calculated  values 
of  s  are  1.21,  1.31  and  1.41  and  1.39  at  x  =  0.45,  0.85,  0.135  and  1.95 
m  from  the  source  for  a  smooth  boundary.  Figures  5-8d,e,f  indicate  that 
the  maximum  concentrations  at  x  =  2.8  m,  X  =  0.49m,  x  =  3.9  m,  X  =  0.68, 


and  at  x  -  6.05  m,  X  =  1.06,  are  no  longer  at  ground- 1  eve  1 . 


rsi 


The  heights  of  the  maximum  concentration  of  the  vertical  profiles 
downwind  from  the  elevated  sources,  Figure  5-8g,  are  initially  at  the 
height  of  the  source,  hut  rapidly  approach  the  ground.  Figure  5-8h  shows 
for  example,  that  it  is  impossible  to  tell  whether  the  concentration 
profile  at  X  =  0.49  and  zS/h  =  0.133  is  due  to  an  elevated  or  a  ground- 
level  source.  This  pattern  is  also  observed  in  neutral  flow,  however, 
the  vertical  profile  at  X  -  1.06,  Figure  5 - 8 i  ,  shows  that  the  position 
of  the  maximum  concentration  is  beginning  to  rise  above  the  ground. 

The  same  pattern  is  observed  for  the  unstable  flow  above  the  rough 
boundary . 

The  variation  of  the  measured  maximum  ground-level  concentration 

-  2 

C(x,0,0)Uh  /Q  is  plotted  in  Figure  5-9  for  the  smooth  boundary  and  in 
Figure  5-10  for  the  rough  boundary.  As  seen  from  these  figures,  the 
ground-level  concentrations  from  an  elevated  source  exceed  after  a 
certain  distance  the  ground- level  concentrations  from  lower-level 
sources.  Comparison  with  Figure  5-11  shows  a  very  close  similarity 
between  the  calculated  and  the  measured  diffusion  pattern,  and  the 
dimensionless  concentration  values  in  the  two  figures  are  of  the  same 
order  of  magnitude,  although  not  identical.  Comparison  of  the  rough  and 
smooth  floors  data  shows  that  the  increased  shear  has  slightly  decreased 
the  maximum  ground-level  concentration  for  the  source  at  7.  /h  =  0.133 
and  has  slightly  modified  the  values  of  the  dimensionless  concentrations 
at  different  point;,  however,  the  basic  pattern  of  diffusion  has  not 
c  hanged . 

Figures  5-12  and  5-1.3  show  the  cross-wind  spread  of  the  plume  and 


its  dependence*  on  the  height  of  the  source.  It  can  be  seen  that  is 
larger  for  g  round- 1  eve  I  sources,  as  predicted  by  I.amb  (19  79).  The 


measured  values  of  o  / h  are,  however,  larger  than  the  calculated  values 
(see  Lamb,  1981,  Figure  5.3). 

Figure  5-14  compares  the  measured  values  of  the  cross-wind  spread 
over  the  smooth  floor  with  the  Prairie  Grass  data  and  the  results  of 
Deardorff  and  Willis  (1975)  and  Lamb  (1979).  It  appears  that  the  wind- 
tunnel  measurements  at  small  values  of  xw"“/(hU),  give  larger  values  of 
(Jy/h.  However,  they  are  consistent  with  some  Prairie  Grass  measure¬ 
ments.  At  larger  distances  the  wind-tunnel  data  is  in  better  agreement 
with  the  few  available  field  measurements. 

The  difference  between  the  values  of  a  /h  for  short  distances  could 

y 

be  due  to  the  effect  of  the  surface  layer  present  in  the  wind-tunnel 
experiment . 
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Figure  5-1.  Temperature  profiles,  smooth  floor 
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Figure  S-2. 


Mean  velocity  profiles,  smooth  floor. 
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Figure  5-4.  Vertical  distribution  of  the  vertical  velocity  fluctuations  (w')  /(w* 
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Figure  5-8.  Selected  dimensionless  concentration  profiles  C(x,0,z)  Uh  /Q,  for  smooth  floor  (SF) 
zS/h  =  0.133  and  rough  floor  (RF)  zS  =  0  (continued). 
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Figure  5-10.  Dimensionless  ground-level  concentrations,  rough  floor 
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6.  CONCLUSIONS  AND  COMMENTARY 

Measurements  of  mean  and  turbulent  velocities  and  temperatures  over 
smooth  and  rough  floors  for  stable,  neutral  and  unstable  flows  in  the 
MWT  at  CSU  were  found  to  be  similar  to  those  measured  in  many  typical 
boundary  layers  corresponding  to  the  appropriate  stability. 

The  wind-tunnel  stable  boundary  layer  flows  had  a  Monin-Obukhov 
length-scale  of  approximately  16  cm  for  the  smooth  floor  and  25  cm  for 
the  rough  floor.  In  the  unstable  boundary  layer,  the  corresponding 
length  scales  were  7  cm  and  16  cm  respectively. 

Diffusion  patterns  in  the  neutral  boundary  layer  from  a  ground- 
level  source  were  similar  to  those  found  in  the  Marchwood  Engineering 
Laboratory  experiments  by  Robins  (1978),  and  these  measurements  were 
also  in  good  agreement  with  Pasquill's  diffusion  estimates.  Diffusion 
data  in  the  vertical  plane  from  the  elevated  source,  for  both  surface 
roughnesses  considered,  showed  that  near  the  source  the  maximum  con¬ 
centration  occurred  at  the  height  of  the  source,  then  it  rapidly 
approached  the  ground  with  increasing  distance  from  the  source. 

In  stably  stratified  turbulent  flows,  wind-tunnel  data  showed  that 
with  increasing  distance  from  the  source,  vertical  diffusion  of  the  gas 
was  much  less  than  under  neutral  stratification.  This  is  easily  ex¬ 
plained  since  increasing  vertical  displacement  in  stable  flows  requires 
a  steady  input  of  energy,  assuming  that  the  density  of  the  diffusing  gas 
remains  unchanged.  Hence  in  stable  flows,  with  increasing  distance  from 
the  source,  the  plume  center  remains  at  approximately  a  constant  height, 
equal  to  the  release  height  of  the  gas. 

The  patterns  of  diffusion  from  ground-level  and  elevated  sources  in 
unstable  flows,  for  both  the  smooth  floor  and  the  rough  floor,  were 
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similar  to  that  found  in  the  water-tank  experiment  of  Deardorff  and 
Willis  (1975)  and  in  numerical  models.  The  measurements  confirmed  the 
predicted  rise  of  plumes  from  ground-level  sources  after  a  travel  time 
on  the  order  of  0.5  h/w*  and  an  initial  rapid  descent  of  plumes  from 
elevated  sources.  They  also  showed  that  after  a  certain  distance  from 
the  source,  ground-level  concentrations  from  an  elevated  source  were 
higher  than  ground-level  concentrations  from  sources  at  lower 
elevations . 

The  measured  cross-wind  spreads,  a^,  were  found  to  be  in  agreement 
with  the  Prairie  Grass  measurements.  The  measurements  confirmed  the 
numerical  predictions,  that  is  larger  for  ground-level  sources  than 
for  elevated  sources. 

Thus,  it  may  be  concluded  that  the  wind-tunnel  measurements  confirm 
the  predicted  unique  features  of  diffusion  in  a  CBL,  which  cannot  be 
described  by  a  simple  Gaussian  model  or  by  a  conventional  Eulerian  type 
mass  balance  differential  equation.  They  also  indicate  the  CBL  diffu¬ 
sion  pattern  can  be  found  for  values  of  -h/L  on  the  order  of  6  which  are 
often  encountered  in  the  atmosphere.  The  effect  of  surface  roughness  or 
the  diffusion  pattern  is  found  to  be  small. 

The  work  exhibits  the  value  of  physical  simulations  in  wind 
tunnels,  which  can  be  used  to  model  stratified  flows  with  much  more 
complicated  boundary  conditions  as  well  as  diffusion  of  buoyant  and 
negatively  buoyant  flows  in  ASL.  Such  problems  cannot  easily  be 
simulated  numerically  (Lamb,  1981). 
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cross-sectional  area  of  source  L 

d 

diameter  of  particles  L 
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buoyancy  flux  L^T  ^ 

g 
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acceleration  due  to  gravity  LT 

h 

height  of  source  L 

a  ^ 

plume  rise  L 

•  Iz 

constant  of  Equation  (3) 

SL 

scale  of  turbulent  eddies  L 

A  Q 

3  -1 

volumetric  discharge  rate  L  T 
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relative  fall  velocity  of  particles  =  Vg/U 

r 

o 

radius  of  plume  at  source  L 

■ 

characteristic  response  length  of  particle  L 

diffusion  Lime  T 

u 

characteristic  mean  velocity  of  air  LT  * 

M  0 

mean  velocity  of  air  at  effective  stack  height  LT  * 

V 

vo 1 ume  L 

vg 

particle  fall  velocity  LT 

!•  « 

exit  velocity  of  plume  LT 

•  x,y,z 

a  right-hand  coordinate  system  with  origin  at  the  L 

base  of  the  simulated  point  source 

X 

•  • 

z 

0 

position  of  maximum  deposition  rate  of  particles  L 

aerodynamic  surface  roughness  length  L 

p 

ratio  of  Lagrangian/Eulerian  time  scale 

p 

9 

-3 

mass  density  of  plume  ML 

• 

>■  ; 
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tV  •  •  .  •  .  •  •  •  -  v  .  v-  .  .  ■ 

‘-1  *-J»- .1-0  *- 

L 


model  length  scale 


model  velocity  scale 
7^  standard  deviation  of  the  vertical  plume  spread  L 

r  particle  relaxation  time  T 


L  -  length 
r  -  time 
ll  -  mass 
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WIND-TUNNEL  SIMULATION  OF  PARTICLE  PLUMES  IN  THE 
ATMOSPHERIC  SURFACE  LAYER 

One  of  the  objectives  of  the  first  year's  program  was  to  explore 
possible  methods  of  modelling  the  dispersion  of  stable  aerocol loidal 
particle-plumes  in  the  wind  tunnel  and  to  design  (or  the  basis  of  this 
investigation)  the  second  year's  research  program  (Tasks  2-A  and  2-B) . 

The  results  of  this  investigation  are  presented  herein. 

1 .  Classification  of  Heavier-than-Air  Particle  Plumes 

The  following  analysis  focuses,  initially,  on  steady-state,  stable, 
particle-plumes  (PP)  from  ground- level ,  or  elevated  sources.  Two  param¬ 
eters  which  significantly  affect  the  behavior  of  such  plumes  are:  the 

relative  fall  velocity  of  the  individual  particles  r  =  V  /U,  where  V 

8  g 

is  the  particle  fall  velocity  and  U  the  characteristic  mean  velocity 
of  the  air;  and  the  relative  negative  buoyancy  of  the  air-particle 
mixture . 

Several  types  of  particle-plumes  which  can  be  classified,  according  to 
these  two  parameters  follow. 

•  Particle-plumes  with  both  negligible  relative  particle  fall 
velocity  and  negligible  plume  buoyancy.  The  dispersion  of 
such  PP  is  identical  to  that  of  passive  tracer  plumes  and,  as 
shown  iii  previous  investigations  (Cermak,  1981  and  Snyder, 
1981),  they  can  be  satisfactorily  simulated  in  meteorological 
wind  tunnels. 

Local  concentrations  of  the  tracer  gas  in  such 
simulations  are  proportional  to  the  spatial  probability 
density  function  of  particles  dispersing  from  the  same  source, 
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and  local  flux  values  cu  *  ds  are  proportional  to  the 
probability  of  particles  passing  through  the  area  ds. 

This  limiting  case  will  be  designated  as  the  Passive 
Reference  Plume  (PRP).  Its  dispersion  characteristics  at 
different  categories  of  atmospheric  stability  and  surface 
roughness  were  simulated  in  Task  1-B  of  the  project. 

When  the  concentration  of  particles  is  increased,  the  plume 
may  be  exposed  to  buoyancy  effects,  even  though  the  relative 
fall  velocity  of  the  individual  particles  remains  negligible. 
(This  study  refers  only  to  negatively  buoyant  plumes,  although 
it  is  theoretically  possible  to  obtain  buoyant  plumes  using 
lighter-than-ai.’.  particles.)  Such  plumes  behave  like  nega¬ 
tively  buoyant  gas  plumes  and,  as  shown  in  previous  investiga¬ 
tions  by  Snyder  (1981)  and  Poreh  (1981),  they  can  be  fully 
simulated  in  meteorological  wind  tunnels  by  keeping  the  dimen¬ 
sionless  buoyancy  flux  constant  in  the  model  and  atmosphere. 

This  type  of  plume  is  designated  as  a  Buoyant  Plume  (BP). 
Its  dispersion  characteristics  are  scheduled  for  investigation 
during  the  second  year  of  the  project  (Task  2-A) . 

It  is  not  easy  to  produce  large  buoyancy  effects  in  the 
atmosphere  with  large  concentrations  of  particles.  However, 
they  may  be  produced  from  a  mixture  of  particles  and  a  carrier 
gas  whose  density  is  less  than  that  of  air. 

A  different  limiting  case  is  obtained  when  the  concentration 
of  the  particles  in  the  plume  remains  small.  However,  the 
relative  fall  velocity  of  the  individual  particles  is  no 
longer  negligible.  Each  particle  will  experience  a  gravity 


force  which  causes  a  mean  downward  motion  and  an  early  deposit 
on  the  ground.  No  interaction  between  particles  is  expected 
in  this  case  and  motion  of  the  particles  is  not  expected  to 
modify  the  atmospheric  flow.  Thus,  the  time  averaged  plume 
characteristics  are  also  equal  to  the  ensemble  average  of  a 
single  particle  dispersion.  This  type  of  plume  will  be 
referred  to  as  a  Hejivy  Particle  Plume  (HPP)  .  Only  monodis- 
persed  HPP  should  be  investigated  and  characterized ,  whereas 
the  PRP  and  the  BP  could  be  composed  of  different  size  par¬ 
ticles,  so  long  as  the  relative  fall  velocity  of  the  largest 
size  particles  remains  negligible. 

•  When  only  a  small  fraction  of  a  polydi spersed  particle-plume 
has  a  non-negl igib le  relative  fall  velocity,  that  fraction 
will  disperse  as  a  HPP,  whereas  the  rest  of  the  particles  will 
disperse  as  a  PRP.  When  the  rest  of  the  plume  behaves  as  a 
BP,  modifications  to  the  velocity  field  may  occur,  partic¬ 
ularly  near  the  source,  and  the  dispersion  of  that  fraction, 
which  is  affected  by  the  velocity  field,  will  differ  from  that 
of  the  HPP.  This  is  called  an  HPP/BP  combination. 

Figure  A- 1  provides  a  graphic  relationship  of  the  four  defined 

particle-plume  regimes  in  terms  of  relative  fall  velocity  (V  /U)  and 

8 

3  6 

dimensionless  buoyancy  flux  (  F/(xu'  )  x  10  ). 

2.  Required  Relative1  Part  ic  U-  Fall  Velocities 

The  properties  of  an  HPP  depend  on  the  relative  particle  fall 
velocity,  r.  When  r  is  small,  the  buoyancy  effect,  which  distin¬ 
guishes  the  HPP  from  the  PRP,  will  not  be  detectable.  Thus  it  is  of 
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primary  interest  to  estimate  the  relative  fail  velocities  at  which  the 
effect  can  be  detected. 

figure  A-2  shows  the  effect  of  a  0.4  downward  shift  of  a 

Gaussian  plume  on  the  vertical  concentration  distribution.  It  appears 
from  the  figure  that  positive  identification  of  the  differences  between 
HPP  and  PRP  characteristics  require  that  the  downward  shift,  D,  satify 
t  he  i nequa 1  i  ty 

D  >  0.4  a  .  ( A- 1 ) 

z 

The  downward  shift  of  the  plume  is  expected  to  be  of  the  order  of 

D  ~  V  •  t  =  V  x/U,  (A-2) 

g  g 

where  t  is  the  diffusion  time  and  U  is  the  average  velocity  of  the 
plume,  taken  as  the  velocity  at  the  height  of  the  source. 

In  order  to  satisfy  Eq.  2  it  is  required  that 

r  =  V  /U  >  0.4  a  /x  =  0.4  I  .  (A-3) 

g  -  z  z 

The  ratio  I  is  a  mild  function  of  x.  However,  for  our  estimates,  an 

z 

average  constant  value  of  may  be  used  for  each  stability  category. 

These  values  were  estimated  using  modified  estimates  of  o ^  (Pasquill, 
1974,  p.  375)  for  zq  =  10  cm  at  x^  =  1  km  (approximately  1.66  m  in  the 
wind  tunnel)  and  are  shown  in  Table  A-l,  together  with  the  calculated 
values  of  r  .  This  table  also  contains  the  corresponding  fall  veloci¬ 
ties  and  diameters  of  unit  density  spheres  descending  in  still  air  at 
those  fall  velocities,  assuming  U  =  2.0  m/ s . 

The  downward  motion  of  the  HPP  affects  the  distribution  of  the 
ground-level  concentrations.  Stewart  (1967)  compared  his  experimental 
data  to  estimates  obtained  from  the  "Classical  Theory"  and  "Statistical 
Theory"  methods  for  predicting  the  diffusion  of  particles.  The 
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"Classical  Theory"  appears  to  predict  more  accurately  the  position  of 

the  maximum  deposition  rate  of  particles,  x  ,  whereas  the  "Statistical 

m 

Theory"  provies  better  estimates  of  the  deposition  rate. 

Figure  A-3,  from  Stewart  (1967),  shows  the  theoretical  estimates  of 

x  /h,  where  h  is  the  height  of  the  source,  for  different  values  of  I 
m  z 

and  r  .  The  figure  indicates  that  the  value  of  r  =  0.03  would  give  a 

10-15  percent  decrease  of  x  / h  for  unstable  flows  (B)  .  A  similar 

m 

change  would  occur  for  r  =  0.015  and  0.01  in  neutral  and  stable  flows. 
These  results  confirm  the  previous  estimates  of  the  minimum  relative 
particle  fall  velocities  required  for  obtaining  detectable  HPP 
characteristics  in  wind-tunnel  experiments  shown  in  Table  A-l. 


3.  Criteria  for  HPP  Simulations 


Dispersion  in  wind-tunnel  models  requires  that  the  ASL  flow  is 
correctly  simulated.  It  is  evident  from  the  preceding  discussion  that 


an  additional  requirement  for  the  simulation  of  HPP  is 


V 

_g 

u 


\u  , 


(A-4) 


\  /  ">  \  /  P 

It  is  usually  not  difficult  to  satisfy  this  requirement.  The  ratio  I ^ 
should  also  be  equal  in  the  model  and  the  atmosphere.  This  condition  is 
satisfied  when  the  ASL  is  correctly  simulated. 

To  simulate  the  dispersion  of  heavy  particles,  it  is  also  required 
that  the  particles'  response  to  the  turbulent  eddy  motion  be  correctly 
modelled.  When  a  free  falling  particle  encounters  a  small  change  in  the 
ambient  air  velocity  it  will  adjust  to  the  new  conditions  within  a  time 
of  the  order  of 


T  =  Vs/*  , 


( A-5 ) 
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which  is  designated  as  the  particle  relaxation  time .  During  this  time 
the  particle  will  descend  a  vertical  distance  of 


S  =  V8  ' 


(A-6) 


The  length  scale,  s,  may  be  termed  the  characteristic  response  length  of 
the  particle. 

To  achieve  a  similar  response  of  the  particles  in  the  model  it  is 
required  that  the  dimensionless  number 


v  2/g 


£  A.  V8Vp 


(A-7) 


where  £  is  a  typical  scale  of  the  turbulent  eddies.  The  scale  of 

eddies  in  the  model  is  proportional  to  the  linear  scale  of  the  model  and 

thus  it  follows  that  when  the  model  scale  is  X-  =  £  /£  =  L  /L  ,  all 

T  m  p  m  p 

the  velocities  in  the  model  should  scale  as 


.  _  8m  Um  .  1/2 

'll  V  ~  U  \ 


(A-8) 


This  requirement  implies  a  considerable  reduction  of  the  model  Reynolds 

3/2 

number,  which  is  scaled  down  as  A  .  The  model  Reynolds  number  may 

not  be  reduced  below  a  critical  value  and  thus  the  above  requirements 

can  be  met  only  in  large  models,  or  when  simulating  high  wind  speeds. 

At  high  wind  speeds,  however,  most  particle  plumes  behave  as  PRP  and 

such  cases  are  of  no  interest  in  the  proposed  study. 

2 

Fortunately,  when  /g £  becomes  very  small  (namely  the  response 

length  of  the  particles  is  small  compared  to  the  length  of  the  eddies), 
the  dispersion  of  particles  would  be  approximately  equal  to  that  of 
"tagged"  fluid  particles  and  independent  of  particle  inertia.  Using  the 
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analysis  ot  t'sanady  11961),  as  summarized  by  Pasquill  (1974),  p,  152, 
K  *  | .  3-133,  it  is  estimated  Lh.jt  when 


R 


1 

In 


0.16 


( A-9 ) 


(where  2  is  the  smallest  significant  eddy  size  in  the  flow),  the 
particle's  motion  will  be  independent  of  its  inertia  and  will  fully 
respond  to  the  turbulent  motion. 

Consider  diffusion  from  an  elevated  source.  The  size  of  the  energy 
containing  eddies  is  of  the  jrder  of  h,  where  h  is  the  height  of  the 
source  above  ground.  The  energy  contained  in  eddies  whose  size  is 
smaller  than  0.1  h  is  small  and  their  effect  on  the  diffusion  can  be 
neglected.  Thus,  Eq.  A-9  with  £  =  0.1  h,  or 

V  2 

-Jh  <  0.015  ,  (A- 10) 

is  a  necessary  condition  to  ensure  that  the  particle  inertia  is  not 
affecting  diffusion. 

According  to  Smith's  (1961)  analysis  on  the  growth  of  a  cluster  of 

particles  descending  in  a  turbulent  flow,  the  size  of  the  cluster,  a, 

would  decrease  when  the  fall  velocity  V  became  large,  due  to  the  fact 

that  the  particles  were  crossing  the  turbulent  eddies.  Smith’s  analysis 

suggests  that  for  large  times  (see  Pasquill,  1974,  p.  15,  Eq.  3.1.30), 

the  reduction  in  cluster  size  for  small  V  /U  would  be  of  the  order  of 

8 

2  2  2 

V  (S  / ( 4  IP),  where  is  the  ratio  of  the  Lagrangian/Eulerian  time 

scales  (see  Pasquill  1974,  p.  135,  Eq.  3.89).  When  (3  =  5 

V 

UR  <  0.125  (A- 1 1  ) 


and  the  effect  of  the  fall  velocity  on  a  is  less  than  10  percent. 
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The  effect  of  the  fall  velocity  on  continuous  PP  is  expected  to  be 
smaller  than  its  effect  on  clusters  of  particles,  but  in  view  of  the 
uncertainties  in  Lhe  analysis  and  the  estimate  of  p,  the  requirement  of 
F.q .  A- 1 1  will  be  adopted  in  this  analysis.  It  follows  that  when  both 
Eq.  A-9  and  A-ll  are  satisfied,  the  diffusion  of  particle  plumes  is 
independent  of  inertia  and  response  time  and,  therefore,  Eq.  A-8  need 
not  be  satisfied.  However,  the  limitations  imposed  per  Eqs .  A-9  and 
A-ll  must  be  satisfied  in  both  model  and  atmosphere. 

As  evidenced  in  Figure  A-4,  many  cases  of  atmospheric  HPP  will 
satisfy  Eq.  A-ll  and  usually  also  satisfy  Eq.  A-10,  as  well.  However, 
the  simulated  HPP  will  satisfy  the  latter  equation  only  if  the  scale  of 
the  model  is  not  too  large.  Figure  A-4  reveals  the  dependence  of  the 
scale  reduction  factor,  1/A^,  on  the  relative  fall  velocity,  r,  of  the 
particles  and  the  wind-tunnel  speed,  U,  at  an  elevation  corresponding  to 
a  60  m  high  prototype  source.  The  figure  reveals  an  inverse  relation¬ 
ship  between  r  and  1/A^,  for  a  given  wind-tunnel  velocity.  Maximum  size 
of  the  model  is,  however,  limited  by  the  physical  dimsnions  of  the  wind 
tunnel . 

The  average  ABL  is  typically  600  m  high  (Robins,  1978).  In  the  MWT 
it  is  possible  to  generate  a  shear  layer  1.2  m  thick  so  that  with  a 
scale  of  1:500,  the  entire  ABL  can  be  modelled.  However,  it  is  usually 
sufficient  to  create  a  shear  layer  whose  thickness  is  four  times  the 
height  of  the  source.  Therefore,  a  minimum  scale  of  1:200  will  suffice 
to  model  a  60  m  tall  prototype  source  in  the  MWT.  This  constraint  is 
also  shown  in  Figure  A-4.  In  spite  of  the  multiple  limitations  it 
appears  from  the  figure  that  a  variety  of  HPP  could  satisfactorily  be 
simulated  in  the  Meteorological  Wind  Tunnel. 
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4 .  Negatively  Buoyant  Particle  Plumes 

When  the  combined  density  of  the  air-particle  mixture  leaving  the 
source  is  larger  than  that  of  the  ambient  air,  the  plume  will  descend 
below  the  source  height  in  exactly  the  same  way  as  a  hot  buoyant  plume 
would  rise  above  the  stack  height. 

It  is  important  to  examine  the  conditions  at  which  such  a  descent 
would  be  detectable,  both  in  wind-tunnel  experiments  and  in  the  atmo¬ 
sphere.  The  effect  can  be  estimated  using  Briggs  (1969)  plume  rise 
equation  (see  Pasquill,  1974,  pp.  242  and  248), 

AH  =  1.8  F1/3  x2/3  U_1  (A-12) 

where 

F  =  U)rQ2  g  Ap/p  =  Qg  Ap/(np)  (A- 13) 

is  the  buoyancy  flux  (divided  by  7t) ,  rQ  -  the  radius  of  the  initial 
jet,  u)  -  the  exit  velocity  and  Q  -  the  (volumetric)  discharge  at  the 
origin. 

It  is  necessary  to  determine  for  which  conditions 

^  >  0.4  ,  (A-14) 

z 

at  x  =  1.67  m  and  x  =  1000  m,  in  neutral  flows,  where  a  /x  is 
m  p  z 

approximately  0.04.  It  follows  from  Eq.  A-14  that  the  above  conditions 
require  that 

>  0.4  x  10"6  .  (A- 15) 

Ux 

Using  the  values  U  =  1.5  m/sec  and  x  =  1.67  m,  it  is  found  that 
6  m 

F  >  4.0  x  10  6  m^/sec3  ,  (A-16) 

m 

whereas  for  the  prototype, 
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(A- 17) 


F  >  2.4  x  10  '  in  ’/sec 
P 

By  assuming  a  6  mm  orifice  for  a  source  diameter  and  establishing  w  and 
U  =  1.5  m/sec  the  source  strength,  Q  =  uiA  becomes 


Q  =  42  x  10  ^  m'Vsec  =  0.25  lit/min 


(A-18) 


The  required  density  difference  in  the  wind  tunnel  thus  becomes 


.  F  •  n 

J —  =  o.o3 

p  Qm  •  g 


(A-19) 


The  density  ratio  of  glass/air  is  approximately  2000.  Therefore 


the  relative  volume  fraction  of  the  glass  particles  is 
AV/V  =  15  x  10"6  ’ 


(A-20) 


and  only  0.015  cc,  or  0.033  gram,  of  glass  need  be  added  to  each  liter 


of  air. 


The  linear  void  ratio  of  the  mixture  in  this  case  is 


1  -  (AV/V) 1/3  =  1  -  0.024  S  0.976 


(A-21) 


so  that  the  "separation"  between  particles  will  be  approximately  40 
diameters . 

The  above  buoyancy  flux  was  calculated  for  AH =  0.4,  which  is 
estimated  to  be  the  minimum  descent  which  can  accurately  be  measured  in 
the  wind  tunnel.  It  would  probably  require  that 


*»  <  0.16 
CT 

Z 


(A-22) 


to  completely  eliminate  negative  buoyancy  effects.  Since  AH  is 

1/3 

proportional  to  F  ,  it  would  be  necessary  to  reduce  Ap/p,  or  F,  by 
a  factor  of  15  to  achieve  this  requirement.  Similarly,  to  increase  the 
effect  by  a  factor  of  m,  the  buoyancy  flux  would  have  to  be  increased 


by  m  . 


Appendix  A 


128 


In  most  atmospheric  PP,  the  particles  will  descend  due  to  their  own 
fall  velocity  and  not  by  modifying  the  density  of  the  plume.  Buoyancy 
effects  can  also  be  created  by  a  mixture  of  particles  and  a  heavier,  or 
a  lighter,  carrier  gas. 

5 .  Conclusions 

The  various  types  of  PP  are  classified  according  to  their  relative 
fall  velocities  and  dimensionless  buoyancy  flux: 

PRP  -  Passive  Reference  Plume 

BP  -  Buoyant  Plumes 

HPP  -  Heavy  Particle  Plumes 

HPP/BP  -  Heavy  Particle  Plume  diffusing  together  with  a  Buoyant 

Plume . 

All  four  types  of  plumes  can  be  simulated  in  the  wind  tunnel.  The 
simulation  of  HPP  is  restricted,  however,  to  a  finite  range  of  relative 
fall  velocities  which  depend  on  the  scale  of  the  model  and  the  height  of 
the  source. 

Figure  A-l  schematically  describes  the  regions  where  the  four  types 
of  plumes  occur  for  neutrally  stable  flows. 
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Table  A-l.  Minimum  Settling  Velocities  for  Hl’F 


Stabi lity 

a^/x 

V  /u 

8 

v  (n 

8 

V2 

Category 

(at  1  km) 

(-) 

cm/  sec 

pm 

B 

0.080 

0.032 

6 . 4 

47 

D 

0.038 

0.015 

3-0 

32 

E 

0.023 

0.010 

2.0 

26 

(1)  Assuming  U  =  2.0  m/sec. 

(2)  Assuming  U  =  7.0  m/se..  and  a  unit  density  sphere. 

(3)  Assuming  U  =  2.0  m/sec  and  0.25  gm/cc  density 
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Figure  A-2.  The  Effect  of  a  Vertical  Descent  D  = 
on  the  Concentration  Distribution 
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X 


X 


m 


Figure  A-3.  Comparison  of  Theoretical  Estimates  of  the  Location 
of  Maximum  Deposition  for  Various  Values  of  r  and 
Iz,  from  Stewart  (1967). 
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Scale  Reduction  Factor  -l/X 


Vg/  U 


Figure  A-4.  Constraints  tot  Wind-Tunnel  Si  mu  J  n  t  i  o>-,  (,,,  j  HPP  Source 

Height  60  m.  Neutral  ASi, .  ( I.  x.imp  I  r :  U  -  2.0  m/sec, 

V  / U  -  0.06,  200  <  A  •  600 

g  h 
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LIST  OF  SELECTED  SYMBuLS 

Symbol 

Definition 

Dimens i on 

A(s , n) 

function  of  shape  factor  and  velocity  profile 
exponent,  defined  by  equation  39 

- 

B(s,n) 

function  of  shape  factor  and  velocity  profile 
exponent,  defined  by  equation  41 

- 

C 

Gaussian  concentration  function 

-3 

ML 

C 

max 

maximum  concentration 

ml"3 

CWIC 

cross-wind  integrated  concentration 

ML"2 

f  (z) 

dimensionless  concentration  profile 

- 

h 

source  height 

L 

Im 

mth  moment  of  dimensionless  concentration  profile 

Lm+1 

Jm 

flux  related  moment 

Lm+2T_1 

Ky 

lateral  mass  diffusivity 

l2t_1 

Kz 

vertical  mass  diffusivity 

l2t_1 

n 

exponent  of  velocity  profile 

- 

Q 

source  height 

mt”  1 

s 

shape  parameter 

- 

U* 

shear  velocity 

lt"1 

U 

mean  velocity 

lt"1 

UH 

reference  velocity 

lt"1 

x,y  ,z 

a  right-hand  coordinate  system  with  origin  at  the 
base  of  the  simulated  point  source 

L 

7. 

centroid  of  concentration  profile 

L 

flux  distribution  centroid 

ZR 

reference  length 

L 

« ( s  ,  a ) 

function  of  shape  factor  and  velocity  profile 
exponent,  defined  by  equation  46 

r 
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P  constant  (=  £ n2  =  0.693) 


A  height  at  which  C/C  =  0.5 

z  max 

flux  standard  deviation 

a  standard  deviation  of  horizontal  plume  spread 


ct^  standard  deviation  of  vertical  plume  spread 

a  centroid-centered  standard  deviation 

zc 

a  estimated  standard  deviation 

ze 


L  -  length 
T  -  time 
M  -  mass 
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LENGTH-SCALES  DESCRIBING  THE  VERTICAL 
DIFFUSION  FROM  GROUND-LEVEL  SOURCES 


1 .  Introduction 

The  vertical  growth  of  passive  plumes  diffusing  from  ground-level 
sources  measured  in  wind-tunnel  simulations  appears  to  have  a  consider¬ 
ably  lower  rate  of  growth  than  that  given  in  widely  used  predicting 
schemes,  like  the  Pasquill-Gif ford  (PG)  cr  curves,  as  demonstrated  in 
Figure  B-l  taken  from  Robins  (1978). 

It  will  be  shown  that  part  of  this  difference  may  be  due  to  the  use 
of  different  length-scales  for  describing  the  vertical  size,  ct^,  of  the 
plume.  Length-scales  used  in  the  literature  to  describe  diffusing  plumes 
will  be  reviewed  and  correlated  analytically,  making  it  possible  to 
compare  different  data  sets  and  to  improve  predictions  of  ground-level 
concentrations . 

2 .  Simple  Gaussian  Model  and  the  Sigma  Curves 

A  solution  of  the  diffusion  equation  for  a  homogeneous  field. 


0 


ac 

3x 


+  K 


where  the  velocity  U,  and 
constants,  is  the  Gaussian 


the  eddy  di f fusivi ties 
concentration  function, 


K 


y 


and  K 

z 


(B-l) 


are 


2 


2 


C(x,y,z)  =  — — -  exp(-  _L  )  exp(.  ),  (B-2) 

2na  a  U  2a  2a 

y  z  y  z 

where  h  is  the  height  of  the  source,  Q  is  the  source  strength,  0  and  o 

y  z 

are  the  standard  deviations  of  the  concentration  profile  in  the  y  and  z 
direction.  The  values  of  a  (x)  and  a  (x)  in  this  solution  are 

y  z 


a 

y>z 


( 2K  x/U) ^ 
y  ,z 


(B-3) 


Appendix  B 


143 


where  h  is  the  height  of  the  source. 

The  pioneering  field  experiments  in  the  British  Meteorological 

Office  in  1958  (Pasquill,  1974)  indicated  that  neither  o  nor  follow 

the  x  law  (Eq.  B-3) ,  and  therefore,  empirical  curves  for  predicting  the 

lateral  and  vertical  spreads  of  plumes  were  proposed.  Perhaps  most 

common  are  the  PG  o  and  a  curves  which  were  adopted  in  Turner's 

y  z 

"Workbook  of  Atmospheric  Dispersion  Estimates"  (1969).  Several  improved 
curves  for  different  stability  categories  and  surface  conditions  have 
appeared  since  (Briggs,  1973),  but  the  so-called  Simple  Gaussian  Model 
(SGM)  (Eqs.  B-4,  B-5  and  B-6)  remains  the  basis  for  calculation  of  the 
concentration  field  from  the  a  and  o  curves. 

y  z 

One  of  the  known  limitations  of  the  SGM  is  its  use  of  a  constant 
characteristic  velocity  U  to  describe  the  flow  field.  In  cases  of 
elevated  sources  it  is  customary  to  assume  that  0  equals  the  mean 
velocity  at  the  effective  height  of  the  source, 
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or  the  mean  velocity  through  the  vertical  extent  of  the  plume,  (h-2 o  ) 
<  z  <  (h+20^)  (Turner,  1969),  which  is  usually  close  to  u(h). 

No  accepted  convention  has  been  adopted  for  ground-level  sources. 
Turner  (1969)  remarks  that  "the  surface  wind  is  most  applicable  to 
surface  or  low-level  emissions,  especially  under  stable  conditions". 
Indeed,  the  most  widely  used  value  is  probably  that  of  u  at  z=10  m, 
which  usually  appears  in  meteorological  data.  However,  the  value  of  u 
at  z  =  2  m,  and  other  heights,  is  also  used.  Several  authors  choose  a 
characteristic  velocity  at  the  "effective  height"  of  the  plume  in  the 
particular  experiment  or  application.  In  practice  only  one  value  of  U 
is  used  in  the  calculations,  even  when  this  recommendation  is  followed 
(see  for  example  problems  1  and  2  in  Turner's  workbook). 

Inherent  problems  make  it  difficult  to  measure  the  vertical 
distribution  of  the  concentration  in  field  experiments  and  thus  the 
vertical  plume  standard  deviation  has  usually  been  estimated  from 
Eq.  B-5,  namely: 


CT 

ze 


Q 

no  C  U 


(B-8) 


y  max 

The  length-scale  obtained  in  this  procedure  is  denoted  by  a^,  rather 

than  o  ,  and  will  be  referred  to  as  the  Estimated  Standard  Deviation  to 
z  - - 

distinguish  it  from  the  o  value  obtained  by  direct  calculation  of  the 
second  moment  of  the  actual  vertical  concentration  profile. 

Different  values  of  0  will  obviously  alter  the  value  of  o  . 
Moreover,  since  Eqs.  B-l  through  B-5  are  mass-balance  equations  for  a 
uniform  flow  field,  the  use  of  a  constant  value  of  0,  which  does  not 
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vary  with  x,  violates  the  mass-balance  equations  and  large  errors  could 

occur  if  a  a  (x)  curve  calculated  from  field  data  at  short  distances 
ze 

were  extrapolated  to  longer  distances. 

3 .  Vertical  Distribution  of  Concentrations 

Wind-tunnel  simulation  of  turbulent  diffusion  from  ground-level 
line  sources  by  Shlien  and  Corrsin  (1976)  revealed  a  deviation  in  the 
vertical  concentration  profile  from  the  Gaussian  model,  as  did  similar 
studies  by  Poreh  and  Cermak  (1964).  Malhotra  and  Cermak  (1964)  also 
found  similar  deviations  in  wind-tunnel  simulations  of  diffusion  from 
point  sources  in  both  neutral  and  unstable  conditions.  Their  data  fit 
the  form, 

C(z)  =  C  exp(-  p(z/A  )S],  (B-9) 

Ilia  X  Z 

where  p  =  £n(2)  =  0.693  and  A  is  the  height  where  C/C  =  0.5,  with  a 
value  of  s  =  1.4.  Similar  experiments  in  neutral  flows  by  Robins  (1978) 
fit  the  same  form  with  s  =  1.7.  Field  observations  (Pasquill,  1974,  pp. 
204-205)  also  indicated  that  Eq.  B-9,  with  s  <  2,  described  the  measured 
vertical  profile  better  than  the  Gaussian  curve  (s  =  2).  The  Porton 
experiment  gave  s  =  1.15  at  x  =  100  m,  whereas  the  Cardington  experi¬ 
ments  gave  s  =  1.5  at  x  =  229  m.  Both  experiments  were  conducted  in  a 
neutral  atmosphere.  The  Prairie  Grass  experiments  suggest  that  s  =  1.48 
±  .28  at  x  =  100  m.  The  higher  values  of  s  correlate  with  fairly 
stable  atmospheric  conditions. 

Several  solutions  of  the  diffusion  equation  with  variable  u(z)  and 
K^(z),  reviewed  by  Pasquill  (1974)  and  Robins  (1978),  also  show 
departures  from  the  Gaussian  form. 

Thus,  it  appears  that  the  concentration  field  downstream  from  a 
point  source  can  be  described  in  the  form, 
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C  =  C  exp(-  y2/2CT  2)  f(z),  (B-10) 

max  r  3  y  ’ 

or  the  equivalent  form 

C  =  Cmax  exp['  P(y/Xy)2j  f(z)’  (B-ll) 

where 

f(z)  =  exp[-  P(z/\z)Sl,  (B-12) 

with  s  varying  in  the  range  1  <  s  <  2. 


4 .  Vertical  Plume  Length-Scales 

Various  length-scales  are  used  to  characterize  the  vertical  size  of 
plumes.  These  are  defined  using  the  dimensionless  concentration  profile 
f(z)  =  C(z)/C  and  its  moments  I  ,  where 

ID3X  m 

00 

I  =  /  zm  f(z)  dz  .  (B— 13) 

m  0 

The  one-half  peak  concentration  height,  is  defined  by 

f (A  )  =  0.5  .  (B-14) 

z 

Similarly,  heights  corresponding  to  other  fractions  of  the  peak  have 
been  used,  such  as  the  one-tenth  peak  h  ..  ht.  The  one-half  peak  height 
is  easier  to  determine  because  of  the  large  gradient  at  that  height. 

The  standard  deviation,  ct^,  is  the  most  commonly  used  vertical 
measure  of  the  plume.  It  is  defined  by 

o2  =  (yi/  .  (B-15) 

The  centroid  of  the  concentration  profile ,  z,  is  defined  as 

z  =  I1/IQ  .  (B— 16) 
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i'he  centroid-centered  standard  deviation,  a  ,  is  defined  as 
-  -  - - >  zc 


■\2 


Ozc  =  S  (z-z)  f(z)  dz  /IQ  , 


(B-17 


and  it  is  related  to  a by 


2  2-2 
a  =  a  -  z  . 
zc  z 


(B-18 


The  length-scales  z  and  o  ,  according  to  Shlien  and  Corrsin,  are 


important  because  z  is  approximately  equal  to  the  mean  particle  dis¬ 
placement  perpendicular  to  the  wall,  and  the  second  centroid-centered 
moment  is  approximately  equal  to  the  dispersion  of  the  particles  around 
the  centroid. 


The  previously  defined  Estimated  Standard  Deviation,  a  ,  was 
r  J - ♦  ze> 


calculated  from  the  ground-level  concentration  (Eq.  B-8).  A  similar 
length-scale  is  defined  by 


h  = 


CWIC  •  u*  ’ 


(B-19) 


where  u*  is  the  shear  velocity  and  CWIC  is  the  Cross-Wind  Integrated 
Concentration,  defined  as 


CWIC  =  /  C(y,0)  dy  . 
-00 


(B-20) 


Since  C(y)  is  well  described  by  a  Gaussian  form,  it  follows  that 


and 


CWIC  =  J2h  C  a 
max  y 


(B-21 


h  =  Jl/h  (U/u*)  aze  .  (B-22 

The  advantage  of  using  h,  rather  than  Oz>  is  that  the  shear  velocity 
is  clearly  defined.  However,  it  is  more  difficult  to  determine  shear 
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velocity  than  it  is  to  determine  the  mean  velocity,  a<_  an  agr.  .d  upon 
height,  say  10  m. 

Using  the  vertical  concentration  profile  given  by  Eq.  B- 11,  the 

four  length-scales  a  ,  z,  a  and  A  can  be  correlated. 

7.  zc  7. 


00 


I 

m 


m 


2 


expl-£(z/Az)s] 


dz 


=  p- (m+1)/s  a  m+1  s'1  r 

Z 


[ (m+l)/s] 


(B-23) 


and  thus, 


Y  =  P'1/s  ir(3/s)/r(i/s)]^ 


(B-24) 


and 


Y  =  P'1/s  r(2/s)/r(i/s)  . 

z 

The  value  of  a  /A  can  be  calculated  from 
zc  z 


<  r5  >2  -  <  x-  > 


2 


(B-25) 


(B-26) 


z  z  z 

Table  B-l  and  Figure  B-2  show  the  dependence  of  the  above  ratios  on 
the  shape  parameter,  s. 

5 .  Flux  Related  Length-Scales 

The  length-scales  z(x)  and  ^(x)  are  related  to  the  concentration 
distribution  at  the  distance  x  and  thus  give  information  on  the 
"residing'’  matter  at  that  position.  The  moments 


00 

J  =  J  zm  u(z)  f(z)  dz 

m  rv 


(B-27) 


on  the  other  hand,  are  related  to  the  distribution  of  the  flux  at  each 
distance,  normal  to  the  mean  wind-direction,  and  may  be  useful  in 
describing  the  motion  of  clouds  and  particles. 
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The  moment  .]  is  proportional  to  the  total  flux,  Q 


J  =  J  u(z)  f(z)  dz  =  Q/(V27T  o  ). 
0  y 


(B-28) 


Hie  flux  distribution  centroid  is  given  by 


Zf  =  Ji/Jo 


(B-29) 


and  the  flux  standard  deviation  is  given  by 


°c  =  (JJJ, 


(B-30) 


f.  2  0  ' 

To  calculate  the  above  moments  and  length-scales,  the  vertical  velocity 
distribution,  u(z),  must  be  known.  Two  forms  of  describing  this 
distribution  are  widely  used:  the  Law  of  the  Wall  representation,  which 
has  a  sound  theoretical  basis;  and  the  power-law  representation,  which 
is  often  used  in  practical  applications  since  it  is  more  convenient  for 
calculations.  Estimates  of  the  power  n  for  various  surface  and 
stability  classes  were  made  by  Irwin  (1979)  and  are  shown  in  Table  B-2, 
taken  from  Hanna  (1982). 

Using  the  power  law 


y~  =  (  ?_  )n 

U  K  z  J  ’ 

R  R 


(B-31) 


where  is  a  convenient  reference  height  within  the  region  described  by 
the  power  law,  the  various  flux  related  moments  become 


rn  n 
ZR 


00  U 

f  m+n  c \  ,  R 
j  z  f(z)  dz  =  — 

0  zj 


m+n 


(B-32) 


L  fo 1  I ows  that 


■(l+n)/s  1+n 


“  =  “  P 


C  , J2n  o  z  ' 

max  y  R 


s"1  T  [ ( l+n)/s  1  (B-33) 
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=  P  [r  (  ^  )  /  r  (  )] 


(B-34) 


^  =  P  s  (  r  (  ^5  )  /  r  (  ^5  )i*^ 

Z 


(B-35) 


_f  =  r[(3+n)  /  s]  r  (1/s) 
az  r(3/s)  r  [ (l+n)/s ] 


(B-36) 


1  =  r  K2+n)/s]  r  (i/s) 

:  r  (2/s)  r  [(i+n)/sj  * 


(B-37) 


The  dependence  of  these  ratios  on  s  and  n  is  also  shown  in 
Table  B-l  and  Figure  B-2.  Using  Eqs.  B-24  and  B-33,  the  flux,  Q,  can  be 
written  as 


Q  s  n  \  \  Ur  Snax  *  <VV  A(s’n)  ’ 


(B-38) 


where 


where 


A(s,n)  =  0.9493  yf2/n  p“(1+n)/s  s"1  r[(l+n)/s]  ,  (B-39) 


Q  =  71  °y  CTz  UR  Cmax  B(s’n)  * 


B(s,n)  =  +2,n  s~  r  -LLLtOl^.L 

ir(3/s)  /  rci/s) ) 1 1  n)U 


(B-40) 


(B-41) 


The  dependence  ot  A(s,n)  and  B(s,n)  on  s  and  n  is  shown  in 
Table  B-l  and  Figures  B-3  and  B-4,  respectively. 

The  relationship  between  and  o  can  calculated  directly 

from  Eqs.  B-8  and  B-40, 


a  U  a 

ze  _  R  .  z  .n  ,  . 

_  _  (  )  B(s ,n) 

z  U  ZR 


( B— 42 ) 
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Since  Uz~  is  a  constant  in  a  boundary  layer  satisfying 


Eq.  B-31 , 


o  o 

ze  /  z  \  n  /  \ 

—  =  1  r.  >  B(s’n)  - 

z  U 


(B-43) 


where  z-  is  the  height  where  U  is  being  measured,  or 


a  z  7—-  -1/ (1+n)  z  n/(ln)  -l/(l+n) 

-5-  =  (  )1+n  B  =  -S  B 

a  ct  ct 

ze  ze  ze 


(B-44) 


The  above  relationship  clearly  reveals  that  for  a  given  s,  n  and 

/v  -  —  — »  -  1  4-  —  /r  l/(l  )  nn  f  kiwi  f  k  n  »/-'>♦  .  /vf  44  4/1  ri  f  U  /.  4 


z- ,  a  is  proportional  to  o 
U  z  ze 


and  thus  the  rate  of  growth  of 


with  x  is  expected  to  be  smaller  than  that  of  aze>  which  is 
compatible  with  the  measurements  shown  in  Figure  B-l. 

Equation  B-43  may  also  be  written  in  the  form 


o  a 

ze  _  ,  z  .  n 

ct  '  a  -  3 

Z  Zy 


(B-45) 


where 


a(s,n)  =  B(s,n) 


(B-46) 


Equation  B-45  may  be  used  to  determine  the  height,  relative  to  o ■  , 
at  which  the  effective  velocity,  U,  should  be  measured  to  ensure  that 
CTze  =  ct^.  Unfortunately,  or(s,n)  varies  considerably  with  both  s  and 
n.  For  example,  a(1.5,  0.15)  =  4.42,  a(1.5,  0.20)  =  6.12  and  a(2.0, 

0.2)  =  1.69.  Thus,  it  appears  that  measuring  U  at  CTz,  or  at  a 
constant  fraction  of  it,  will  not  ensure  the  above  desired  equality. 


6 .  Improving  Concentration  Predictions 

Given  semi-empirical  values  for  °z(z)  for  different  stability 
categories  and  surface  roughnesses,  it  is  essential  to  determine  whether 
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these  va  I  ut  s  rep  resent  t  hr  true  standard  devia 


>f  the  plume 


concentration  or  estimated  values  obtained  o  nd-level  concent  ra¬ 
tions  and  an  eftritive  velocity  measured  as  i<  height. 

Ground- 1  eve  1  toment  rat  urns  can  be  cal<  directly  from  the 

estimated  values  of  u  using  F.q.  B-5.  Elevated  concentration  levels, 
particularly  for  z  ■  o  should  be  estimated,  however,  using  Ec| .  R-lf 
To  carry  out  such  a  procedure  it  is  necessary  to  estimate  both  n  and  s. 

Given  estimates  of  o  ,  Eu .  B-40  should  be  used  to  determine  G 

7.  max 

at  ground- 1  eve  1 ,  using  the  velocity  at  any  reference  height,  except  very 
close  to  the  ground  where  the  power  law  is  not  expected  to  hold  It  is 

necessary,  of  course,  to  estimate  both  s  and  n.  The  value  of  n  can 
be  estimated  using  Table  B-2.  No  procedure  tor  estimating  s  is  avail¬ 
able,  but  even  a  rough  estimate  of  the  shape  parameter  will  usually  give 
a  more  accurate  prediction  than  one  obtained  by  the  simple  Gaussian 
model . 
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Table  B-l.  Ratios  of  Different  Length-Scales  as  a  Function  of  s  and  n 


n 

O 

i/A 

o  /A 
zc 

O 

/A 

* 

r/A 

0 

c/K 

O 

0  fo 
zc 

0 

r'* 

0,c/0f 

V°2 

z 

t'z 

0.  /o 
i  c  zc 

A 

B 

.05 

1 

20 

1 

473 

1 

085 

.996 

1 

522 

1 

134 

1 

U15 

1 

357 

,67b 

1 

.  6u7 

1.033 

1 

044 

] 

019 

.848 

.665 

.05 

1 

40 

1 

189 

0 

902 

.775 

1 

225 

0 

939 

0 

7  86 

1 

318 

.651 

1 

304 

642 

1.030 

1 

04  1 

1 

015 

.781 

.766 

05 

1 

60 

1 

025 

0 

794 

.648 

1 

054 

0 

825 

0 

656 

1 

290 

.  b32 

1 

278 

.622 

1.028 

1 

038 

1 

012 

.740 

848 

.05 

1 

80 

0 

921 

0 

725 

.567 

0 

944 

0 

751 

0 

573 

1 

269 

.616 

1 

257 

.606 

1.026 

1 

036 

1 

010 

.713 

T915 

.05 

2 

00 

0 

849 

0 

678 

.512 

0 

870 

0 

701 

0 

516 

1 

253 

603 

] 

24 1 

.  593 

1.025 

1 

034 

1 

008 

.694 

.9-70 

.05 

2 

20 

0 

798 

0 

644 

.472 

0 

817 

0 

665 

0 

475 

1 

241 

.592 

1 

229 

.582 

1.024 

1 

033 

1 

006 

.681 

1.016 

.  10 

1 

20 

1 

473 

1 

085 

.996 

1 

570 

1 

182 

1 

033 

, 

157 

.676 

1 

328 

.658 

1.066 

1 

089 

1 

038 

.834 

.641 

.  10 

1 

40 

1 

189 

0 

902 

.775 

1 

260 

0 

975 

0 

798 

1 

*18 

.651 

1 

292 

.633 

1.059 

1 

081 

1 

030 

.762 

.742 

.  10 

1 

60 

1 

025 

0 

794 

.648 

1 

082 

0 

854 

0 

663 

1 

290 

.632 

1 

266 

.613 

1.055 

1 

075 

1 

024 

.718 

.823 

.  10 

1 

80 

0 

921 

0 

725 

.567 

0 

968 

0 

777 

0 

578 

1 

269 

.616 

1 

246 

.597 

1.052 

1 

071 

1 

019 

.690 

.889 

.  10 

2 

00 

0 

849 

0 

678 

.512 

0 

891 

0 

724 

0 

520 

1 

253 

.603 

1 

231 

.583 

1.049 

1 

068 

1 

015 

.670 

.944 

.  10 

2 

20 

0 

798 

0 

644 

.472 

0 

836 

0 

685 

0 

478 

1 

241 

.592 

1 

219 

.572 

1.047 

1 

065 

1 

012 

.656 

.989 

.  15 

1 

20 

1 

473 

1 

085 

.996 

1 

617 

1 

229 

1 

051 

1 

357 

.676 

1 

316 

.650 

1.098 

1 

133 

1 

056 

.823 

.621 

.  15 

1 

40 

1 

189 

0 

902 

.775 

1 

295 

1 

Oil 

0 

809 

1 

318 

.651 

1 

280 

.625 

1.089 

1 

121 

1 

044 

.746 

.720 

.  15 

1 

60 

1 

025 

0 

794 

.648 

1 

109 

0 

863 

0 

670 

1 

290 

.632 

1 

255 

.604 

1.082 

1 

112 

1 

035 

.700 

.801 

.  15 

1 

80 

0 

921 

0 

725 

.567 

0 

991 

0 

802 

0 

583 

1 

269 

.616 

1 

236 

.588 

1.076 

1 

105 

1 

.027 

.669 

.867 

.  15 

2 

00 

0 

849 

0 

678 

.512 

0 

911 

0 

746 

0 

523 

1 

253 

.603 

1 

221 

.574 

1.072 

1 

100 

1 

021 

.648 

.921 

.  15 

2 

20 

0 

798 

0 

644 

.472 

0 

853 

0 

70S 

0 

480 

1 

241 

.592 

1 

210 

.563 

1.069 

1 

.096 

1 

017 

.633 

.966 

.20 

1 

20 

1 

473 

1 

085 

.996 

1 

665 

1 

277 

1 

068 

1 

357 

.676 

1 

304 

.642 

1.130 

1 

.176 

1 

073 

.815 

.603 

.20 

I 

40 

1 

189 

0 

902 

.775 

1 

329 

1 

046 

0 

819 

1 

318 

.651 

1 

270 

.616 

1.117 

1 

160 

1 

057 

.733 

.701 

.20 

1 

60 

1 

025 

0 

794 

.648 

l 

136 

0 

912 

0 

677 

1 

290 

.632 

1 

245 

.596 

1 . 108 

1 

148 

1 

04S 

.683 

.781 

.20 

1 

80 

0 

921 

0 

725 

.567 

1 

014 

0 

826 

0 

58? 

1 

269 

.616 

1 

227 

.579 

1 . 101 

1 

139 

1 

036 

.651 

.846 

.20 

2 

00 

0 

849 

0 

678 

.512 

0 

930 

0 

767 

0 

526 

1 

253 

.603 

1 

212 

.566 

1.095 

1 

132 

1 

026 

.629 

.900 

.20 

2 

20 

0 

798 

0 

644 

.472 

0 

871 

0 

725 

0 

483 

l 

241 

.592 

1 

201 

.554 

1.091 

1 

.127 

1 

021 

.613 

.946 

.25 

1 

20 

1 

473 

1 

085 

.996 

1 

712 

1 

324 

1 

085 

1 

357 

.676 

1 

293 

.634 

1.162 

1 

.219 

1 

090 

.809 

.587 

.25 

1 

40 

1 

189 

0 

902 

.775 

1 

363 

1 

081 

0 

829 

1 

.318 

.651 

1 

260 

.608 

1.146 

1 

.198 

1 

070 

.721 

.684 

.25 

1 

60 

1 

025 

0 

794 

.648 

1 

163 

0 

940 

0 

683 

1 

290 

.632 

1 

236 

.5  88 

1.134 

1 

184 

1 

055 

.669 

.  763 

.25 

1 

80 

0 

921 

0 

725 

.567 

1 

036 

0 

850 

0 

591 

1 

269 

.  b!6 

1 

218 

.571 

1.125 

1 

172 

1 

043 

.635 

.829 

.25 

i 

00 

0 

849 

0 

678 

.512 

0 

950 

0 

789 

0 

529 

1 

253 

.603 

1 

204 

.557 

1.118 

1 

164 

1 

033 

.611 

.883 

.25 

2 

20 

0 

798 

0 

644 

.472 

0 

888 

0 

744 

0 

484 

1 

241 

.592 

1 

193 

.545 

1.113 

1 

157 

l 

025 

.594 

.928 

.30 

1 

20 

1 

473 

1 

085 

.996 

1 

758 

1 

370 

1 

102 

1 

357 

.676 

1 

283 

.627 

1.194 

1 

263 

1 

106 

.805 

.573 

.30 

1 

40 

1 

189 

0 

902 

.775 

1 

396 

1 

116 

0 

839 

1 

318 

.651 

1 

251 

.601 

1.174 

1 

237 

1 

083 

.712 

.669 

.30 

1 

60 

1 

025 

0 

794 

.648 

1 

169 

0 

968 

0 

690 

1 

290 

.632 

1 

228 

.580 

1.160 

1 

219 

1 

064 

.656 

.748 

.30 

1 

80 

0 

921 

0 

725 

.567 

1 

057 

0 

874 

0 

596 

1 

269 

-bl6 

1 

210 

.563 

1.149 

1 

205 

1 

050 

.620 

.813 

.30 

2 

00 

0 

849 

0 

678 

.512 

0 

968 

0 

809 

0 

532 

1 

253 

.  603 

1 

197 

.549 

1.140 

1 

194 

1 

039 

.595 

.867 

.30 

2 

20 

0 

798 

0 

644 

.472 

0 

905 

0 

763 

0 

486 

1 

241 

-592 

1 

186 

.537 

1.133 

1 

186 

1 

029 

.578 

.912 

.40 

1 

20 

1 

473 

1 

085 

.996 

1 

850 

1 

463 

1 

133 

1 

357 

.676 

1 

265 

.612 

1.256 

1 

348 

1 

138 

.803 

.550 

.40 

1 

40 

1 

189 

0 

902 

.775 

1 

462 

1 

184 

0 

857 

1 

318 

.651 

1 

235 

.586 

1.229 

1 

312 

1 

106 

.698 

.645 

.40 

1 

60 

1 

025 

0 

794 

.648 

1 

240 

1 

023 

0 

701 

1 

290 

.632 

1 

212 

.565 

1.209 

1 

287 

1 

082 

.636 

.723 

.40 

I 

80 

0 

921 

0 

725 

.567 

1 

100 

0 

920 

0 

603 

1 

269 

.616 

1 

196 

.548 

1.195 

1 

268 

1 

063 

.596 

.788 

.40 

-* 

00 

0 

849 

0 

678 

.512 

1 

005 

0 

850 

0 

537 

1 

253 

.603 

1 

183 

.534 

1.183 

1 

254 

1 

048 

.568 

.84] 

.40 

2 

20 

0 

798 

0 

644 

.472 

0 

537 

0 

799 

0 

469 

1 

241 

.592 

1 

172 

.522 

1.174 

1 

242 

1 

036 

.549 

.886 

.50 

1 

20 

l 

473 

1 

085 

.996 

1 

941 

1 

555 

1 

163 

1 

357 

.676 

1 

249 

.599 

1.318 

1 

432 

1 

168 

.809 

.533 

.50 

1 

40 

1 

189 

0 

902 

.775 

1 

526 

1 

251 

0 

874 

1 

318 

.651 

1 

220 

.573 

1.283 

1 

38 

l 

129 

.691 

.627 

.50 

1 

60 

1 

025 

c 

794 

.648 

1 

290 

1 

0,6 

0 

712 

1 

290 

.632 

1 

199 

552 

1.258 

1 

354 

I 

098 

.621 

.  705 

50 

1 

80 

0 

921 

0 

725 

.567 

1 

141 

0 

964 

0 

610 

1 

269 

.616 

1 

183 

.534 

1.239 

1 

330 

1 

075 

.577 

.769 

.50 

2 

00 

0 

849 

0 

678 

.512 

1 

040 

0 

886 

0 

541 

1 

253 

.603 

1 

171 

.520 

1.225 

1 

311 

1 

057 

.547 

.822 

.  50 

2 

20 

0 

798 

0 

644 

.472 

0 

969 

0 

834 

0 

492 

1 

241 

592 

1 

161 

.508 

1.213 

1 

296 

1 

041 

.525 

.867 

60 

1 

20 

I 

473 

J 

086 

.996 

2 

031 

1 

64S 

1 

m 

1 

357 

.6  76 

i 

235 

.587 

1.379 

I 

516 

1 

297 

.822 

.521 

.60 

l 

40 

1 

189 

0 

902 

.  775 

1 

589 

1 

316 

0 

890 

1 

318 

.651 

1 

207 

.560 

1.336 

1 

459 

1 

149 

.68ft 

.614 

60 

I 

hO 

1 

025 

0 

704 

,64ft 

1 

338 

1 

127 

0 

721 

1 

290 

.632 

1 

187 

.539 

1.306 

1 

419 

1 

113 

.62  l 

.691 

.60 

i 

80 

0 

921 

0 

725 

.56  7 

l 

161 

1 

007 

0 

61b 

1 

269 

.  b  16 

1 

172 

.522 

1.283 

1 

389 

1 

Gfib 

.5  62 

.755 

bO 

2 

00 

0 

849 

0 

678 

.512 

1 

074 

0 

926 

0 

5«5 

1 

253 

.  b03 

1 

lbO 

.507 

1.265 

l 

366 

1 

.064 

.529 

.809 

.60 

2 

20 

0 

798 

0 

644 

.472 

0 

999 

0 

868 

0 

494 

1 

241 

.592 

1 

151 

.495 

1.251 

1 

348 

1 

046 

.505 

853 

Table  8-2.  Estimates  of  the  Power  n  in  Eq.  31*  by  Irwin 
from  Hanna  (1982) 


z  (m) 

o 

Class: 

A 

B 

C 

D 

E 

F 

0.01 

0.05 

0.06 

0.06 

0.12 

0.34 

0.53 

0.10 

0.08 

0.09 

0.1 1 

0. 16 

0.32 

0.54 

1.00 

0.17 

0.17 

0.20 

0.27 

0.38 

0.61 

3.00 

0.27 

0.28 

0.31 

0.37 

0.47 

0.69 

*  2  <  100  m. 
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